RNA Biology

Escherichia coli YicC enzyme is the founding member of a family of endoribonucleases that is encoded in virtually all bacterial species. Previous structural studies revealed that this ribonuclease binds RNA by a novel mechanism in which the hexameric apoprotein presents an open channel that undergoes a large rotation upon RNA binding and clamps down on the RNA. The current study follows up on these findings by examining the cleavage of various oligonucleotide substrates designed to probe recognition elements required for YicC binding and cleavage. A 26-nucleotide RNA oligomer (oligo), with a KD in the low micromolar range, was the standard to which numerous oligos with altered sequence were compared. In vitro RNase assays and fluorescence anisotropy binding measurements indicated that the preferred substrates for YicC were relatively small RNAs that contain some secondary structure. Larger RNAs or highly structured RNAs were less-than-optimal substrates. Similarly, RyhB RNA, a [~]90-nucleotide, iron-responsive RNA of E. coli, which has been described as a target of YicC binding and/or cleavage, was a poor YicC substrate in our assays. These results suggest that the native substrates for YicC-family members are very small RNAs or RNA fragments derived from larger RNAs.

Transfer RNAs (tRNAs) are known for delivering amino acids to the growing polypeptide chain during translation. They can also influence gene expression, especially in times of nutrient starvation, through differential tRNA expression and modification. tRNAs have a highly consistent cloverleaf structure, but relatively few known regulatory elements govern this conserved structure despite the 20 different standard isotypes. This study examines gene enrichment patterns near tRNA in 1154 fungal genomes. Genes enriched in proteasome regulation, ion transport, and rRNA were found to be significantly closer to tRNAs than other pathways. These results were consistent across KEGG over-representation analysis (ORA), KEGG Gene Set Enrichment Analysis (GSEA), and Gene Ontology (GO) analysis. Proteasome, ion transport, and RNA are all important aspects of protein production and regulation, suggesting that genes required for the synthesis and quality control of proteins, including tRNAs, are located near each other. Protein regulation is an energetically expensive process, and local co-regulation could increase efficiency and stress impacts on proteins.

Pseudouridine is the most abundant RNA base modification due to its prevalence in tRNA and rRNA, where it serves as a key modulator of structure and function. Yet even in a widely used model organism, the budding yeast Saccharomyces cerevisiae, the positions of all pseudouridines in tRNA have not been completely annotated. Using Nanopore direct RNA sequencing (DRS), an established method for detecting RNA pseudouridylation positions, we sequenced cytosolic tRNA from eight pseudouridine synthase (PUS) knockout S. cerevisiae strains, including deletion strains of Pus1, Pus3, and Pus7. Analysis of these data verified thirty-four existing pseudouridine annotations and predicted eleven previously unannotated pseudouridine sites. Our analysis revealed DRS signal changes at several non-uridine sites with the loss of a PUS, including apparent changes in modification abundances at position 37 upon deletion of Pus3. LC-MS/MS and primer extension assays, however, indicated no change in the abundance of these modifications with the loss of Pus3. Our analysis underscores the need for caution in interpreting DRS-based signal changes, particularly in modification-dense regions. Combining existing modification annotations for the thirty-one isoacceptors in the Modomics database with our dataset that added annotations for the remaining eleven isoacceptors, we created a map of all detected pseudouridines, and the enzymes responsible for their catalysis, across the forty-two S. cerevisiae cytosolic tRNA isoacceptors.

Pseudouridine ({Psi}) is an important post-transcriptional modification of many noncoding RNAs that is under-characterized in microRNA (miRNA) due to historical limitations in pseudouridine mapping methods. {Psi} modification stabilizes RNA duplex structures and could therefore play an important role in miRNA target binding and repression. To investigate the extent to which mammalian miRNAs are modified with {Psi}, we profiled the modification landscape of short (<30 nt) RNA in human cells and mouse tissues using bisulfite sequencing. Our approach was powered to detect small RNA pseudouridylation based on robust detection of known {Psi} positions in tRNA fragments (tRFs), some of which show tissue-specific patterns of modification. In contrast with tRFs, we find that miRNA pseudouridylation is exceedingly rare, with a single modified miRNA (miR-3068-5p) identified in mouse tissues. Pseudouridylated miR-3068-5p diSerentially repressed predicted miRNA targets with less stable miRNA:mRNA pairing modes. This study fills a long-standing gap in transcriptome-wide {Psi} profiling and reveals a new potential function for {Psi} as a modulator of activity of small regulatory RNAs.

Regulation of RNA pools allows for adaptation to changing environments and stress, which is especially important in pathogenic bacteria such as Mycobacterium tuberculosis. RNA degradation is a significant contributor to RNA abundance, and Ribonuclease (RNase) E has a rate-limiting role in degradation of a majority of mycobacterial transcripts. However, many open questions remain about the RNA substrate requirements and specificities for efficient cleavage by mycobacterial RNase E. Here, using both Mycolicibacterium smegmatis and M. tuberculosis RNase E, we demonstrate that this enzyme is only active on substrates with a minimum length of approximately 27 nt. Furthermore, we show that mycobacterial RNase E prefers substrates with 5 monophosphates rather than 5 triphosphates, and that the positions of cleavage events within substrates are dictated by both sequence and distance from the RNA ends. Our results also suggest that RNase E may be affected by product inhibition. Finally, we show that M. smegmatis RNase E behaves similarly to M. tuberculosis RNase E, validating the use of this model organism for RNA degradation studies.

BackgroundSmall non-coding RNAs (sncRNAs) play central roles in post-transcriptional gene regulation. In addition to canonical microRNAs (miRNAs), fragments derived from vault RNAs (vtRNAs), called small vault RNAs (svtRNAs), have been reported in human cells. However, the absence of a standardized annotation framework has hindered their systematic detection, quantification, and comparison across small RNA sequencing (small RNA-seq) studies. MethodsWe developed an expression-based annotation strategy to identify svtRNAs from human small RNA-seq datasets. Using FlaiMapper followed by structure and expression-based filtering, we generated two annotation sets: a stringent "miRNA-like" set enriched in Argonaute-associated datasets, and (ii) a broader "Total" set derived from total small RNA-seq libraries under relaxed structural constraints. We explored the expression of the annotated svtRNAs across the different datasets analyzed: multiple normal and tumor-derived human cell lines, including Argonaute immunoprecipitation datasets. ResultsWe identified a repertoire of svtRNAs that are detected across independent datasets and, in several cases, reach abundance levels comparable to canonical miRNAs. Several highly abundant svtRNAs correspond to molecules with experimental validation from prior studies, supporting the robustness of our annotation strategy. Importantly, the same "dominant" (in terms of gene expression) svtRNAs emerged independently from Argonaute-associated and total small RNA datasets, supporting the idea of enzymatically consistent, reproducible svtRNA processing. We further identified svtRNAs derived from distinct vtRNA precursors that could share identical seed sequences, suggesting the possibility of svtRNA families with potential miRNA-like regulatory properties. We provide a standardized annotation that enables reproducible svtRNA quantification. ConclusionsOur study establishes a comprehensive expression-based annotation resource for human svtRNAs. By enabling their systematic detection and reproducible quantification, we show that svtRNAs appear to represent an abundant component of the human small RNA landscape.

Even though teleost fish and mammals share the same mitochondrial gene content and organization, the teleost mitochondrial transcriptome is still poorly understood. We characterized the mitochondrial transcriptome during zebrafish (Danio rerio) early development by long-read direct RNA sequencing. All heavy-strand specific mRNAs were found to carry 3 poly-A tails of approximately 50-60 residues, and the transcriptome profile was distinctive but practically invariant between stages. Three unusual transcripts were however noted. These included two mRNAs (COI and ND5 mRNAs), with significant 3 untranslated regions corresponding to antisense gene sequences, and a previously not described noncoding RNA named here lncOriL. The ND5 mRNA was found to carry one third of all detected m6A methylation sites in the zebrafish mitochondrial transcriptome. The 313 nt-long lncOriL transcript had an abundance comparable to that of ND5 mRNA and it mapped to mitochondrial genome region covering the origin of light strand replication and four flanking antisense tRNAs. A mitochondrial tRNA-derived fragment (tiRNA5-Asn), with a 35 nt perfect pairing-potential to lncOriL, was present at all stages. Additional analyses including adult zebrafish, scissortail (Rasbora rasbora), and monkfish (Lophius piscatorius) strongly corroborate the results of COI mRNA, ND5 mRNA, and lncOriL transcript prevalence among teleost fish. Surprisingly, our findings in zebrafish were further supported by mitochondrial transcriptome analyses in domestic pig (Sus scrofa) and human (Homo sapiens), including tiRNA5-Asn commonly present in human tissues, suggesting that lncOriL is ubiquitously expressed and regulated in vertebrates. Author SummaryMitochondria contain their own genome and produce essential RNAs needed for energy production. Although fish and mammals share the same mitochondrial gene organization, less is known about how mitochondrial RNAs are processed and regulated in teleost. Using Nanopore direct RNA sequencing, we examined mitochondrial RNAs during early zebrafish development and discovered three unusual transcripts that include extended non-coding regions. Two of these molecules, COI and ND5 mRNAs, carry long 3' untranslated regions formed by antisense gene sequences, suggesting previously unrecognized regulatory potential. We also identified lncOriL, a highly structured long noncoding RNA that spans the origin of light-strand replication and is abundant during development. Strikingly, the same RNA feature, including lncOriL and a matching tRNA-derived small RNA (tiRNA5-Asn), was found not only in zebrafish but also in human mitochondrial transcriptomes. These findings support conservation of regulatory mitochondrial RNAs across main groups of vertebrate species. Our work reveals a new layer of mitochondrial RNA regulation and expands the current understanding of how mitochondrial gene expression is controlled.

Over the past decade, groundbreaking discoveries have cemented transfer RNAs (tRNAs) as versatile regulators of translation and cellular function. As tRNA research gains momentum, several high-throughput sequencing methods for quantitative analysis of tRNA isoacceptors in cells have emerged. However, the strong secondary structure and rich post-transcriptional modification of most tRNA molecules pose significant challenges for reverse transcriptases, thus hampering library preparation and introducing quantification biases. Current approaches rely on processive next generation reverse transcriptases (ngRTs), such as Induro (NEB) and uMRT (RNAConnect), to overcome these problems. Nevertheless, using these commercial enzymes comes with a relatively high cost per reaction. Here, we introduce a redesigned MarathonRT construct with added C-terminal chitin binding domain (CBD) and present a simple and robust one-step purification protocol for the in-house production of the redesigned and the original MarathonRT construct. Next, we set up an affordable colorimetry-based method for determining the specific activity of the enzyme. Our simplified method eliminates protein precipitation and yields over 26,000 enzymatic reactions per 0.5 L culture. Importantly, the in-house produced enzymes showed equal performance to established ngRTs Induro and uMRT in tRNA-seq workflow. In addition, we benchmark the use of rapid tRNA spin column-based extraction method to traditional gel-extraction using tRNA-seq and LC-MS. This improved workflow reduces the time and cost of tRNA-seq library preparation while providing an accessible MarathonRT purification protocol.

Eukaryotes use several distinct quality control pathways to resolve aberrant ribosomes and mRNAs. For example, the no-go decay mRNA pathway is stimulated after ribosome collisions caused by stalled ribosomes translating damaged or truncated mRNAs. Separate decay pathways for non-functional 40S and 60S subunits containing rRNA mutations affecting decoding and peptidyl transferase activity, respectively, have also been elucidated. To our knowledge, whether eukaryotes have evolved a quality control pathway to sense and process globally stalled ribosomes is unclear; however, such a pathway would be advantageous to eukaryotes during exposure to natural elongation inhibitors such as ricin and diphtheria toxin. Here, we test how prolonged robust inhibition of elongation using a high dose of cycloheximide (CHX) affects ribosome turnover. Despite no decrease in cell viability and that mammalian ribosomes have been classically characterized of having a half-life of 3-5 days, a single 24 hr high dose of CHX resulted in drastically shortened half-lives (<24 hr) of 28S and 18S rRNA in A549 cells. A [~]2-fold reduction in nearly all ribosome species was observed by polysome analysis in HeLa and A549 cells after prolonged CHX treatment. Depletion of ribosomes was also evident when assessing ribosomal proteins from both the 40S and 60S subunits by Western blot. Literature supports that ribosomes can be degraded by autophagy and the ubiquitin (Ub)-proteasome system. Upon testing inhibitors of both pathways, only proteasome inhibitors (i.e., MG132 and bortezomib) rescued both rRNA and ribosomal protein levels. Proteasome inhibitors also rescued ribosome levels in polysome profiling experiments. Remarkably, rRNA levels were not rescued during CHX treatment when co-treated with the Ub activating enzyme E1 inhibitor, TAK243. Polysome analysis also showed that the high prolonged dose of CHX did not cause robust accumulation of collided ribosomes compared to control treatments. Proteasome-dependent turnover of rRNA was also observed with high doses of other elongation inhibitors, namely anisomycin, homoharringtonine, and lactimidomycin. The recognition capabilities of the pathway were further expanded as we observed that 80S ribosomes not trapped on the mRNA were also targeted for degradation by the proteasome. Together, our findings define the framework of a regulatory pathway in mammalian cells that degrades both ribosomal subunits in response to prolonged periods of robust elongation inhibition.

RNase R is a processive 3 to 5 exoribonuclease that degrades a broad array of linear RNA species while preserving RNA lariats and circular RNAs (circRNA). In recent years, this enzyme has become pivotal for the field of circRNA research, serving as a key step for circRNA enrichment, purification, and identification. Despite this growing importance, the effects of mutations and truncations in RNase R have been incompletely studied. We make several point mutations and assay their effects on the ability of RNase R to bind and/or degrade RNA substrates. Our data show that selected active site mutations have varying effects on RNA binding and degradation. Furthermore, the increasing interest in circRNA-based RNA therapeutic platforms highlights an urgent need for RNase R in RNA molecular biology labs. However, the substantial cost of commercial RNase R remains a bottleneck, particularly for large-scale studies or the development of circRNA-based technologies. In this protocol, we offer a solution to that problem, namely a more accessible and cost-effective means of purifying high-quality and low-cost RNase R. We provide a highly detailed yet simplified, high-yield protocol that produces recombinant RNase R from Escherichia coli. The method uses a single-step Ni-NTA affinity chromatography procedure without proteolytic tag removal and is optimized for entry-level FPLC systems such as the AKTA Start, ensuring that high-purity enzyme production does not require specialized, high-end instrumentation. A second key feature is the establishment of an optimized reaction framework, including specific buffer compositions and defined enzyme-to-substrate ratios for the purified RNase R. The protocol achieves functional equivalence to premium commercial RNase R, ensuring complete linear RNA digestion without compromising the integrity of circRNA. The combination of a simplified purification workflow and a robust reaction protocol provides an accessible, cost-effective, and reliable solution for any molecular biology laboratory requiring high volumes of RNase R. Key FeaturesO_LIRNase R mutations can block RNase activity, RNA binding or both C_LIO_LIThis protocol purifies [~]40 mg of active RNase R per liter of E. coli culture C_LIO_LIThe protocol avoids medium and high end FPLC systems C_LIO_LIRNase R expression constructs (WT and mutants) will be available on Addgene C_LIO_LIThe protocol includes an optimized reaction buffer to pair with this RNase R C_LIO_LIOptional endotoxin removal step is also included C_LI

The prevalent transcriptional repressor NrdR binds to highly conserved prokaryotic sequences in the promoter regions of operons encoding the essential enzyme ribonucleotide reductase. The NrdR binding sites consist of two partially palindromic 16 bp sequences (NrdR boxes) separated by a 15-16 bp linker sequence. We have assessed the requirement of both boxes for binding, the propensity of different NrdRs to bind to heterologous binding sites, and that the linker sequence is only limited to length and not sequence conservation. As we have observed several deviations from the conserved sequences of the NrdR boxes, we here test the conservation requirements of individual basepairs in the NrdR boxes using a synthetic DNA fragment (Synt DNA) to which the NrdR proteins from the actinomycete Streptomyces coelicolor and the gammaproteobacterium Escherichia coli bind equally well as to their homologous binding sites. By introducing isolated mutations to Synt DNA and testing the binding capacity of NrdR from S. coelicolor and E. coli we expand our understanding of what criteria are needed to build a functional binding site for the NrdR repressor.

BackgroundMetazoan adenosine-to-inosine (A-to-I) mRNA editing temporospatially diversifies the neuronal transcriptome and proteome. The limited read length from next-generation sequencing (NGS) constrains the quantification of the potentially differential editing levels across different splicing isoforms, restricting our understanding of the extent to which RNA editing contributes to molecular diversity and its interplay with splicing. MethodsWe employed reverse transcription nested PCR (RT-nPCR) and developed a novel interfering-Primer PCR (iPrimer PCR) technique to distinguish different transcripts of any gene. We selected multiple essential genes exhibiting RNA editing in coding sequences (CDSs) or untranslated regions (UTRs) for isoform-specific amplification and Sanger sequencing. ResultsNine different Adar isoforms together with pre-mRNA had distinct editing levels at the S>G auto-recoding site, which was predicted to have isoform-specific effects on catalytic activities. Although pre-mRNA editing might exert isoform-dependent promotion/suppression of splicing, closely located editing sites, such as those in neuronal genes qvr and stj, still exhibited high correlation in editing levels due to co-editing. iPrimer strategy further discovered differential recoding levels between the long/short 3UTR isoforms of gene jef. ConclusionsWe provide the first comprehensive solution for isoform-specific PCR amplification of any gene, enabling quantification of RNA editing level of different isoforms. Our results offer insights into how RNA editing interplays with splicing, and highlight its complicated role in expanding molecular diversity. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/725286v1_ufig1.gif" ALT="Figure 1"> View larger version (17K): org.highwire.dtl.DTLVardef@1ebc82org.highwire.dtl.DTLVardef@1ea365dorg.highwire.dtl.DTLVardef@1971aceorg.highwire.dtl.DTLVardef@160d053_HPS_FORMAT_FIGEXP M_FIG C_FIG We developed isoform-specific PCR followed by Sanger sequencing, and achieved the quantification of differential RNA editing levels in different transcripts of a gene.

Accurate detection of RNA splice variants is often hindered when transcripts lack large distinguishable exonic regions, making conventional PCR strategies challenging. We developed a simple melting temperature (Tm)-guided exon-exon junction (EEJ) RT-PCR method to enable variant-specific detection under these conditions. Uni-directional primers spanning exon-exon junctions were designed so that approximately each half anneals to adjacent exons. The Tm of each half-site was set >7{degrees}C below the annealing temperature, preventing stable binding to individual exons and enforcing junction-dependent amplification. The method was evaluated using HTRA1-AS1 long noncoding RNA variants that share overlapping exon sequences but differ in splice connectivity. HTRA1-AS1 comprises five variants, only one with a large distinguishable exon. Tm-guided EEJ primers robustly discriminated the remaining four variants. After optimization, amplification yielded sharp, single bands with minimal cross-reactivity. Compared with conventional designs, this approach reduced heteroduplex and heteroquadruplex formation, improving band clarity. Sanger sequencing confirmed junction specificity, and the method performed well in multiplex settings. Overall, Tm-guided EEJ RT-PCR is a cost-effective, high-resolution approach for detecting RNA variants lacking easily distinguishable exonic regions, readily compatible with standard RT-PCR and qPCR workflows.

Since the discovery of viral Internal Ribosome Entry Sites (IRESes), researchers have sought to find similar elements in mammalian host genes, termed "cellular IRESes". However, the plasmid systems used to measure cellular IRES activity are vulnerable to false positives due to promoter activity in candidate IRESes. Orthogonal methods are needed to validate putative IRESes while carefully avoiding artifacts known to cause false positives. Recently, Koch et al. proposed approaches for studying IRESes, primarily circular RNA-generating plasmids, and for validating mRNA transcripts using smFISH and qRT-PCR. Here, we demonstrate confounding variables and artifacts in each of these approaches that can lead to inappropriate conclusions about potential cellular IRES activity. We show the back-splicing circRNA plasmid creates linear mRNA artifacts associated with false-positive IRES signals. Using orthogonal, gold-standard assays validated with viral IRESes, we find putative cellular IRESes reported using the back-splicing plasmid have no IRES activity. Furthermore, we demonstrate that smFISH and qRT-PCR can misidentify nuclear non-coding RNAs as mRNAs and we validate a single molecule sequencing assay for identifying genuine mRNA 5 ends. Our work establishes reliable methods for robust transcript annotation and IRES studies that avoid documented artifacts arising from bicistronic and back-splicing circRNA plasmid reporters.

Human mitochondrial genome (mtDNA) encodes multiple proteins in the oxidative phosphorylation complexes as well as the ribosomal and transfer RNAs (tRNAs) needed for in situ translation. These genes are transcribed from only three promoters, producing polycistronic transcripts that are co-transcriptionally cleaved by mitochondrial RNase enzymes to release majority of individual gene products. tRNAs separate many of these genes and are thought to serve as "punctuation" marks that enable RNase recognition, binding, and hydrolysis of the 5' "leader" and 3' "trailer" sequences flanking the tRNA. Mutations in the tRNA genes dominate the mtDNA-linked mitochondrial pathologies; yet a systematic study of the impact of tRNA sequence variation on the RNase-catalyzed processing is lacking. Here, we employed human mitochondrial tRNATyr as a model system to dissect the effect of tRNA variants on the in vitro 5' leader and 3' trailer hydrolysis. We found that nucleotide variations located near the catalytic interfaces - particularly within or near the tRNA acceptor stem - showed the strongest defects in 5' processing and prevented release of the downstream tRNA in a tRNA cluster where multiple tRNAs are transcribed in tandem. This work provides mechanistic insight into how mutations disrupt coordinated mitochondrial tRNA processing and establish a framework for predicting variant effects based on their structural position relative to the processing enzymes.

MicroRNAs are short noncoding RNAs that regulate gene expression and are commonly profiled by small RNA sequencing (miRNA-seq). Despite the widespread use of miRNA-seq, datasets are often analyzed with RNA-seq method such as DESeq2 or edgeR, which do not take into account the specific characteristics of miRNA-seq data. Here, we present a benchmark study of normalization and differential expression approaches using a realistic ground-truth dataset. By mixing mouse RNA of two organs, we generated expression trends while capturing biological and technical variability. Using monotonicity across the dataset and expected fold changes from the mixture design, we assessed normalization and differential expression methods. Normalization benchmarking showed that within-sample scaling, particularly Read Per Million (RPM), best preserved the expected monotonic trends, outperforming cross-sample methods such as TMM, rlog, and VST. These approaches sometimes recovered apparent monotonicity among abundant miRNAs, but inspection of individual profiles suggested likely over-correction. Regarding differential expression, edgeR consistently ranked among the best-performing methods across several metrics, including log2 fold-change estimation, with performance comparable to miRNA-seq-specific tools such as miRglmm and NBSR. DESeq2, edgeR-v4, and limma-based approaches tended to systematically underestimate log2 fold changes. Applying a common RPM-based normalization substantially improved the performance of cross-sample methods, highlighting the strong influence of normalization on differential expression analysis. Overall, our findings support within-sample scaling methods such as RPM for normalization, and edgeR, miRglmm, or NBSR for differential expression. The dataset has been made publicly available, providing a valuable resource for objective method comparison and future miRNA-seq software development.

Variants affecting RNA splicing are a major contributor to human disease, yet the consequences of variants outside of the canonical splice motifs are often difficult to determine. Here, we present a protocol for minigene-based evaluation of candidate splice-altering variants. The methodology described includes locus-specific insert design, commercial gene fragment synthesis, and long-read sequencing. The combined approach enables rapid assay development and nucleotide level resolution of the effect on splice isoforms in vitro, providing a scalable framework for functional validation of predicted cryptic splice variants. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=197 SRC="FIGDIR/small/723105v1_ufig1.gif" ALT="Figure 1"> View larger version (42K): org.highwire.dtl.DTLVardef@1a88cb5org.highwire.dtl.DTLVardef@adda98org.highwire.dtl.DTLVardef@1ea587corg.highwire.dtl.DTLVardef@574a63_HPS_FORMAT_FIGEXP M_FIG C_FIG

Muscleblind-like (MBNL) RNA-binding proteins (RBPs) possess modular domains that mediate regulation of alternative splicing and RNA localization. Myotonic Dystrophy Type 1 is a CTG repeat expansion disorder where MBNL is sequestered into intranuclear RNA foci, impairing its function. Previous studies found that MBNL self-associates through its exon 7, but the nature of this interaction is not well understood. We identified a cysteine in MBNL1 exon 7 that enables dimerization through formation of an intermolecular disulfide bond. We likewise demonstrate that MBNL2 dimerizes by forming disulfide bonds between multiple cysteines in its carboxy-terminus. Nucleocytoplasmic fractionation revealed a greater proportion of MBNL1 dimer in the nucleus, suggesting a nuclear function for the MBNL1 dimer. We investigated a connection between MBNL1 dimerization and MBNL1-mediated regulation of alternative splicing. To accomplish this, we mutated the MBNL1 cysteine in question to alanine (C325A) and performed RNAseq. We uncovered novel splicing events sensitive to MBNL1 dimerization. We also found that MBNL1 C325A, when co-expressed with expanded CTG repeats, produces smaller, more numerous foci, suggesting a role for the MBNL1 dimer in maintaining foci integrity. These results provide insight into biological and pathological mechanisms of MBNL1 dimerization and suggest other RBPs might similarly dimerize to regulate function. GRAPHICAL ABSTRACT

Long non-coding RNAs (lncRNAs) are critical regulators of gene expression in eukaryotes. Short reads from Illumina sequencing, reverse transcriptase artefacts, and incomplete second-strand degradation in strand-specific cDNA libraries hamper genome-wide identification of lncRNAs, especially in gene-dense genomes such as Plasmodium. Here, we integrated long-read Oxford Nanopore Technology direct RNA sequencing, ribosome profiling, and single-cell transcriptomics to generate a robust and stage-specific characterization of P. falciparum lncRNAs. We generated comprehensive annotations of lncRNAs expressed in both asexual and sexual blood stages and confirmed their non-coding nature using ribosome profiling. Most lncRNAs showed pronounced stage-specific expression and appeared to be particularly abundant in mature gametocytes. Single-cell RNA sequencing revealed differential expression of many lncRNAs in female and male gametocytes, suggesting important roles in gametocytogenesis and transmission. Many lncRNAs are located antisense to protein-coding genes and are co-expressed with their sense mRNA, possibly from putative bidirectional promoters, while others overlap mRNA coding sequences or 3 untranslated regions and showed negatively correlated expression patterns. Overall, our study shows the prevalence of P. falciparum lncRNAs and highlights their possible roles in controlling the regulation of gene expression, particularly during gametocytogenesis.

Intron length is a fascinating example of form without function. The vast majority of intronic space within genomes remains without a provided utility. It often fascinates us to find introns performing any function at all, establishing an attention bias against the vast lacking of utility of the remaining intergenic space. In an attempt to better understand the greater breadth of intronic length, I investigate here what I term The Kinetic Intron Hypothesis. This hypothesis investigates hypothetical dynamics of intron RNA synthesis and degradation. It explores how NTPs stored within intron RNA might function in mitosis and NTP resource management. Preliminary testing of the hypothesis leads to trends that warrant further exploration and validation by the scientific community. SignificanceCurrently no widely acknowledged model exists to characterize the length of introns within genes, yet intron length is massively abundant in eukaryotic genomes. Here I present an attempt to model the length of introns. In doing so, I explore novel hypothesized intron dynamics, presenting preliminary data for previously uncharacterized intron characteristics. The new data and model have the protentional to unveil new avenues of utility for introns at the intracellular level.